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摘 要 : 运动 想象 脑 电信 号 采集 成 本 高 且 个 体 差异 大 ， 跨 个 体 域 构建 脑 电信 号 模式 识别 模型 属于 典型 的 小 样本 跨 
域 学 习 任 务 。 针 对 该 任务 ， 提 出 了 一 种 运动 想象 脑 电 信号 的 跨 域 特征 学 习 方 法 。 该 方法 首先 选择 最 优 度量 方法 对 
齐 协 方差 并 提取 共同 空间 模式 特征 ; 其 次 ， 在 该 特征 基础 上 采用 领域 自 适 应 方法 学 习 目 标 域 的 最 优 跨 域 特征 。 为 
验证 所 提出 方法 的 可 行 性 与 有 效 性 ， 采 用 经 典 模型 识别 跨 域 特征 ， 在 两 个 公开 的 数据 集 上 进行 对 比 实验 。 实 验 结 
果 表明 ， 通 过 所 提出 方法 学 习 到 的 跨 域 特征 ， 在 运动 想象 模式 识别 中 ， 明 显 优 于 现 有 方法 学 习 到 的 特征 。 此 外 ， 
还 详细 对 比 了 跨 域 特征 学 习 方法 的 各 项 参数 设置 、 性 能 及 效率 。 
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Cross-domain feature learning method for motor imagery EEG signals 
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Abstract: Motor imagery EEG signals requires a high cost for recording, and there is a large difference for individual's signals. 
Cross-subject motor imagery EEG signals recognition task belongs to a typical cross-domain learning problem with a small 
samples set. To solve this problem, this paper proposed a cross-domain feature learning method for motor imagery EEG 
signals to improve the recognition performance. The proposed method first selected the optimal measurement to align the 
covariance of EEG signals, and then extracted common spatial patterns (CSP) from the aligned EEG trials. Second, based on 
the CSP features, an optimal domain adaptation algorithm was selected to learn the optimal cross-domain features for the 
target domain. To validate the feasibility and effectiveness of the learned cross-domain features, a classical model was adopted 
to recognize the learned cross-domain features, and the comparative experiments were conducted on two public datasets. 
Experimental results show that the learned cross-domain features are obviously better than state-of-the-arts methods in 
recognition performance. In addition, this paper also compared the parameters setting, performance and efficiency for the 
proposed method. 
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S 0 ”引言 者 在 没有 肌肉 系统 参与 的 情况 下 控制 轮椅 外 。 近 年 来 ，MI- 
r= x: BCI 已 经 帮助 中 风 患 者 做 康复 训练 以 改善 其 运动 功能 59， 
出 机 接口 BCD02 技 术 绕 过 肌肉 系统 ， 人 允许 大 脑 与 外 部 ”并且 MI-BCI 的 娱乐 性 应 用 已 经 扩展 至 健康 人 群 。 然 而 ， 


x 


设备 之 间 建 立 直接 的 通路 ， 实 现 神经 系统 和 外 部 设备 间 的 直 EEG 信号 记录 人 脑 头 皮层 的 微弱 信号 变化 ， 易 受 噪 声 影响 ， 
接 信息 通信 。BCI 技术 分 为 三 个 主要 步骤 : 收集 大 脑 信号 ， 常用 的 特征 提取 方法 为 共同 空间 模式 (CSP)。CSP 是 应 用 非 
分 析 并 解码 大 脑 信号 ， 再 根据 解码 结果 向 外 部 设备 输出 控制 常 广 泛 的 一 种 空域 滤波 特征 提取 算法 1， 在 运动 想象 任务 
命令 。 头 皮 脑 电 (EEG) 是 BCI 最 常用 的 脑 电 信号 类 型 之 一 ， 的 特征 提取 与 分 类 中 有 不 错 的 效果 ， 该 方法 能 够 从 多 通道 的 
该 信号 采用 非 侵 入 方式 记录 大 脑 活动 的 电 生理 现象 。 基 于 EEG 信和 号 中 提取 出 每 个 类 别 的 最 优 空间 分 布 成 分 。 在 CSP 
EEG 的 BCI 主 要 包括 三 种 方式 : 运动 想象 (MI)， 稳 态 视觉 特征 基础 上 ， 可 采用 线性 判别 分 析 (LDA) 或 支持 向 量 机 
立 (SSVEP)、 事 件 相关 电位 (ERP)。 其 中 ，MI 属于 主动 BCI (SVM) 等 模型 完成 运动 想象 任务 的 模式 识别 。 
模式 ， 因 此 基于 EEG 的 MI-BCI 是 最 常 使 用 的 BCI 范式 ， 实际 上 ， 由 于 EEG 信号 采集 的 成 本 较 高 ， 通 常 单个 被 
也 是 本 文 的 研究 对 象 P]。 试 者 的 EEG 样本 数量 较 少 ， 并 且 每 个 被 试 者 的 EEG 信号 分 
在 MI-BCI 中 ， 用 户 根据 提示 想象 自己 身体 部 位 的 运动 ， 布 存在 个 体 差异 性 ， 无 法 将 跨 被 试 的 全 体 EEG 信号 直接 用 
此 时 会 在 对 应 的 大 脑 皮层 区 域 引起 感知 运动 节律 ， 产 生 大 脑 ” 于 构建 高 性 能 模式 识别 模型 。 因 此 ， 一 些 研究 者 希望 通过 迁 
相关 皮层 区 域 的 脑 节律 调制 ， 从 而 在 EEG 信号 中 编码 不 同 。” 移 学 习 的 方法 降低 EEG 信和 号 特征 分 布 的 个 体 差 异性 外， 并 在 
的 运动 想象 规律 和 特性 。 因 此 ， 不 同 运动 动作 的 想象 可 以 产 。 全 体 跨 被 试 者 的 EEG 信号 上 构建 高 性 能 模式 识别 模型 。 迁 
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寺 
生 不 同 的 空间 局 部 化 解码 ， 从 而 识别 出 不 同 的 运动 感知 节 得 移 学 习 中 最 常见 的 任务 为 领域 自 适 应 ， 经 典 的 领域 自 适 应 方 
模式 ， 最 终 可 根据 模式 识别 结果 的 映射 关系 控制 外 部 设备 。 法 包括 Yang 等 提出 的 迁移 成 分 分 析 (TCA)00、Long 等 提出 
早期 MI-BCI 己 经 应 用 至 医学 领域 ， 用 于 帮助 严重 竣 痪 的 患 ”的 联合 分 布 适 配 方法 JDA)00、 以 及 考虑 稀疏 特征 选择 的 迁 
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录用 定稿 PHAR, €i 运动 


还 包括 Wang 等 提出 的 动 
基 


移 联 合 匹配 改进 方法 (TJMD)D2， 
分 布 自 适应 BDA)03。 在 解决 领域 自 适 应 的 方法 中 ， 其 基 
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Ha iur 由 于 该 方法 需 将 样本 对 齐 至 协 方差 中 
*， 因 此 本 文 将 协 方差 中 心计 算 扩 展 至 不 同 距离 度量 方法 ， 


信号 的 跨 域 特征 学 习 方 法 


目标 是 降低 源 域 和 目标 域 的 特征 分 布 差异 ， 以 保证 源 域 和 
标 域 的 特征 分 布 能 共同 用 于 训练 分 类 器 ， 这 些 方法 可 拓展 
MI-BCI 的 CSP 特征 之 上 ， 用 以 提升 模式 识别 的 性 能 045]。 
在 基于 MI-BCI 的 迁移 学 习 中 ， 一 些 研 究 者 09 结 合 核 共 
同 空间 模式 (KCSP)00 和 迁移 核 学 习 (TKL)08 提 出 了 迁移 核 


想象 脑 
态 
本 
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并 EOF 同 个 体 选 择 最 优 的 协 方差 中 心 度量 方法 。 
1.1 协 方差 矩阵 中 心 对 齐 方法 
协 方差 矩阵 中 心 可 作为 多 维 时 间 序 列 分 布 的 度量 方式 。 
假设 X, 是 被 试 者 进行 运动 想象 任务 获得 的 第 ;个 样本 ， 每 
个 被 试 者 拥有 7 个 EEG 样本 ， 则 个 体 域 的 EEG 信号 样本 协 


共同 空间 模式 方法 (TKCSP)， 并 应 用 于 跨 被 试 者 的 脑 电信 号 
处 理 中 。Jayaraml91 等 研究 者 提出 了 一 种 不 需要 目标 域 标 记 
数据 的 跨 被 试 者 迁移 学 习 的 多 任务 学 习 框 架 ， 并 证 明 该 方法 
对 一 名 肌肉 萎缩 患者 的 迁移 效果 不 错 。AzabL? 等 研究 者 对 
上 述 方 法 做 了 进一步 的 改进 ， 提 出 MI 分 类 中 对 跨 被 试 者 利 
用 Kullback-Leibler 散 度 进行 加 权 迁 移 学 习 。 在 集成 迁移 学 
习 中 ，Hossain[?11 等 研究 者 提出 了 多 类 MI 分 类 中 跨 被 试 迁移 
的 集成 学 习 方 法 。 
在 数据 对 齐 方面 ，ZaniniC25 等 研究 者 提出 了 一 种 黎 曼 距 
离 对 齐 方法 (RA)， 并 采用 改进 的 MDRMP3] 分 类 器 验证 了 该 
方法 的 有 效 性 。YairP4 等 研究 者 提出 了 一 种 在 SPD 和 矩阵 的 锥 
流 形 上 ， 通 过 平行 输 运 解析 表达 式 的 领域 自 适应 方法 ， 并 证 
明 该 方法 在 跨 被 试 者 运动 想象 分 类 、 睡 眠 阶段 分 类 和 心算 识 
别 方面 有 不 错 的 效果 。He 和 WHE RA 的 基础 上 做 了 改进 ， 
提出 了 一 种 欧式 距离 对 齐 方 法 (EA)， 该 方法 在 欧 氏 空间 上 将 不 


方差 矩阵 中 心 碳 为 : -过 xx7 。 该 中 心 采 用 算法 平均 值 表达 ， 
在 经 典 的 EA 算法 中 使 用 欧式 距离 度量 。 随 后 ， 根 据 协 方差 
矩阵 中 心 ， 对 每 个 EEG 样本 进行 对 齐 操作 : 
RY (1) 
经 过 n 次 对 齐 后 ， 所 有 被 试 者 的 EEG 样本 的 平均 协 方 
差 矩 阵 都 等 于 单位 矩阵 ， 因 此 来 自 不 同 被 试 者 的 协 方差 矩阵 
分 布 会 更 为 相似 。 
1.2 ”最 优 协 方差 矩阵 中 心 距 离 度量 方法 
于 不 同 任务 、 不 同 被 试 者 对 EEG 信号 样本 的 分 布 都 
有 影响 ， 在 实际 个 体 域 的 对 齐 中 ， 需 根据 每 个 个 体 域 的 样本 
分 布 情况 ， 灵 活 选 择 合适 的 协 方差 中 心 度量 方法 。 除 了 欧式 
距离 以 外 ， 本 文 将 采用 如 下 7 种 均值 度量 方法 : 
DRZ: 


X -armin Fp, (X. X, ») (2) 


其 中 : d 表示 黎 曼 空间 上 两 个 协 方差 矩阵 的 距离 。 采 用 梯 


同 被 试 者 的 EEG 信号 协 方差 矩阵 对 齐 到 单位 矩阵 上 ， 比 RA 方 
法 更 方便 、 高 效 。 此 外 ，RodriguesPq 等 研究 者 提出 了 黎 曼 


Procrustes 分 析 (RPA) 方 法 ， 用 于 对 齐 跨 被 试 者 的 EEG 信号 ， 该 
方法 在 跨 被 试 MI、ERP 和 SSVEP 分 类 中 有 着 不 错 的 效果 。 
在 流 形 迁 移 学 习 中 ，Zhang 和 Wu27 在 黎 曼 流 形 协 方差 
矩阵 对 齐 基础 上 ， 提 出 了 一 种 在 切 空 间 上 的 流 形 知识 迁移 方 
法 (MEKT)。Singhl”3 等 研究 者 提出 了 一 种 样本 选择 的 迁移 学 
习 方 法 ， 该 方法 从 少量 的 目标 域 样本 中 估计 MDRM 分 类 器 
使 用 的 样本 协 方差 矩阵 。 另 外 ， 基 于 深度 学 习 的 MI-BCI Xf 
移 学习 P9 也 获得 了 不 错 的 效果 。Wmu 等 研究 者 Ba 提出 了 一 种 
应 用 于 MI-BCI 分 类 的 并 行 多 尺度 滤波 器 (CNN)， 提 升 了 分 
类 结果 在 迁移 上 的 性 能 。 
运动 想象 EEG 信号 具有 非 平稳 特性 且 易 受 噪声 影响 ， 
获取 个 体 域 记录 的 信号 代价 高 且 特 征 分 布 差异 大 ， 跨 个 体 域 
构建 脑 电信 号 模式 识别 模型 属于 典型 的 小 样本 跨 域 学 习 任务 。 
从 现 有 研究 的 综述 中 可 以 看 出 ， 目 前 的 研究 集中 在 样本 协 方 


度 下 降 法 ， 在 黎 曼 空间 上 进行 迭代 的 均值 求解 。 


(2) 对 数 欧 式 均 值 : 
X=ep( 1 Znosx | G) 
(3) 对 数 det 均值 : 
Z=- (Xiosx 40.5X, ^ il (4) 
(4) Wasserstein 均值 : 
K (Zaxo) 5) 


其 中 : koX. TRE Wasserstein 度量 方式 ， 初 始 化 迭代 次 
数 ， 不 断 更 新 K ， 直 到 迭代 次 数 停止 


差 对 齐 、 流 形 特征 自 适 应 或 模式 识别 模型 的 迁移 81H。 为 了 
在 个 体 域 之 间 学 习 稳 定 、 可 靠 的 跨 域 特征 ， 本 文 提 出 了 一 种 
运动 想象 脑 电信 号 的 跨 域 特征 学 习 方 法 ， 用 于 提升 跨 个 体 域 
脑 电信 号 模式 识别 性 能 。 首 先 ， 该 方法 选择 最 优 的 协 方差 矩 
车 距离 度量 方法 ， 在 个 体 域 之 间 将 样本 集 协 方差 对 齐 至 单位 
矩阵。 其 次 ， 在 对 齐 后 的 样本 和 集 上 提取 CSP 特征 ， 
JDA、BDA 和 TJM 四 种 不 同 的 领域 自 适 应 方法 ， 在 待 识 别 的 
目标 域 上 学 习 最 佳 的 跨 域 特征 。 本 文 在 两 个 公开 数据 集 上 验证 
了 所 提出 跨 域 特征 学 习 方法 的 可 行 性 与 有 效 性 ， 并 且 从 参数 设 
置 、 模 式 识别 性 能 和 效率 上 进行 了 详细 对 比 实验 分 析 。 

现 有 面向 运动 想象 脑 机 接口 的 跨 域 特征 学 习 方 法 仅 针对 
样本 对 齐 角度 或 CSP 特征 跨 域 自 适 应 角度 ， 跨 域 的 特征 学 
习 能 力 不 够 ， 导 致 运动 想象 模式 识别 性 能 不 足 。 本 文 同时 针 
对 样本 角度 和 CSP 特征 角度 构建 跨 域 特征 学 习 方法 ， 为 个 
体 域 脑 电 信号 样本 和 集 选 择 最 优 的 样本 对 齐 度量 方法 ， 以 及 
CSP 特征 跨 域 自 适应 方法 ， 并 在 两 个 公开 的 脑 电 数 据 集 上 验 
证 了 所 提出 方法 的 模式 识别 性 能 及 效率 。 


1 ”最 优 协 方差 矩阵 样本 对 齐 
在 跨 个 体 域 的 EEG 信号 样本 对 齐 中 ， 本 文选 用 欧式 距 


引入 TCA、 


(5) AJD 对 数 均值 : 
X -Arov[ Xin XB JA (6) 
Wh: MERE B 是 基于 Pham 算法 的 近似 联合 对 角 化 ，4 EB 
的 逆 矩 阵 。 
(6) 调 和 均值 : 
exo) o 


(7) KL 散 度 均值 ， 计 算 方 式 为 
1 n 1 n 1 zi ua 1 n 
Y -(— Wa 12| 二 ES Pais B 
a GA) D» (ix. ) | G2 (8) 


实际 上 ， 该 均值 是 算术 平均 值 和 调和 平均 值 之 间 的 几何 
平均 值 。 
实际 实验 中 ， 本 文 将 根据 不 同 个 体 域 EEG 样本 分 布 情 
况 ， 采 用 trial-and-error 的 方式 选择 最 优 的 度量 方式 ， 计 算 协 
方差 矩阵 中 心算 术 平均 值 ， 并 采用 式 () 对 齐 不 同 个 体 域 的 
EEG 信和 号 样本 集 。 


2 EEG 信号 的 跨 域 CSP 特征 学 习 


在 迁移 学 习 方 法 中 ， 领 域 自 适应 是 常见 的 基于 核 空间 的 
其 本 质 是 最 小 化 源 域 特 


方法 。 征 分 布 和 目标 域 特 征 分 布 的 距 


离 ， 最 常 使 用 衡量 源 域 特征 分 布 和 目标 域 特征 分 布 的 距离 为 
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MMD B21。 在 EEG 信号 的 跨 域 CSP 特征 学 习 中 ， 本 文 的 目标 
是 在 多 个 源 域 和 目标 域 CSP 特征 最 小 化 MMD 距离 。 


3 ”运动 想象 EEG 信号 跨 域 特征 学 习 方 法 


2.1 迁移 成 分 分 析 方 法 3.1 共同 空间 模式 
迁移 成 分 分 析 方 法 (TCA) 为 基于 边缘 分 布 的 自 适应 09 方 k 同 空间 模式 (CSP)B 针 对 二 分 类 任务 ， 从 多 通道 EEG 


法 。 假 设 用 源 域 和 目标 域 的 边缘 分 布 采用 P(X,) 和 P(X,) 表 盲 号 中 提取 可 区 分 特征 。 假 设 X M X; 分 别 为 两 类 不 同 运 动 
示 ，TCA 方法 的 目的 是 减 小 P(X,) 和 P(X,) 之 间 的 差异 。 想象 任务 的 EEG 信号 ， 首 先 ， 计 算 各 域 对 齐 后 EEG 信号 样 


TCA 的 目的 是 最 小 化 mm 距离 ， 其 优化 目标 为 本 集 的 平均 协 方差 矩阵 Ci (i=1,2) ; 其次， 寻找 一 个 最 优 的 
mintr (A7 XLXTA)4 uA 空间 滤波 器 w ,使 得 类 别 1 的 方差 最 大 ， 类 别 2 的 方差 最 小 ， 
v i (9) RRE BP 
SE AAHUA C,w - A(C, C, )w (13) 
其 中 ， 4 是 最 小 化 MMD 的 映射 矩阵， 五 æt ĠO RE, 其 中 : AGB" 分 别 为 广义 特征 值 和 特征 向 量 ， 然 后 利用 


H -L,,,-V (n +m) , 工 即 为 WMWMD 上 距离。 采用 拉 格 朗 日 乘 子 Zun 7 Wu Xy 将 EEG 信号 进行 空域 滤波 获得 特征 向 量 Z o 


法 即 可 求 得 4 的 最 优 解 为 (XLX+y1)"' XHX 最 大 的 严 个 特征 值 对 随后 ， 选 择 特征 向 量 2 的 前 m 行 和 后 m 行 (2m<M)， 形 成 的 

应 的 特征 向 量 。 最 终 ， 通 过 映射 4 变换 源 域 和 目标 域 特征 。 特征 向 量 作为 分 离 后 的 空域 滤波 特征 。 对 于 EEG fü X, 

2.2 联合 分 布 自 适应 CSP 算法 最 终 提取 到 的 特征 f; 79: f -log(ZZT) 。 

联合 分 布 自 适 应 方法 GDA)00 的 目标 是 减 小 源 域 和 目 3.2 线性 判别 分 析 

标 域 的 联合 概率 分 布 的 距离 。 用 P(X,) 和 P(X,) 之 间 的 距离 、 线性 判别 分 析 ILDA)B5 是 MI-BCI 中 常用 的 分 类 器 ， 对 

以 及 P(YIX,) 和 PC|X,) 之 间 的 距离 来 近似 表达 。 由 于 目标 域 提取 到 的 CSP 特征 分 类 后 ， 完 成 运动 想象 任务 的 模式 识别 。 

没有 标记 数据 ， 所 以 无 法 直接 求解 PQ|X,) 。JDA 方法 通过 25 E BEES TE Duyan, Xau 27 分 别 表示 第 ie1{1,2} 
T 迭代 的 方式 在 源 域 特征 和 样本 上 训练 简单 的 分 类 器 ， 并 在 类 运动 想象 EEG 信号 提取 到 的 跨 域 特征 ， 均 值 向 量 和 协 方差 
一 标 域 特征 上 和 预测 伪 标 签 总 。 形 式 地 ，JDA 的 优化 目标 为 矩阵 。 首 先 ，LDA 分 别 计算 为 “类 内 散 度 矩阵 ”$5, 和 “类 间 散 度 
m: 3 HERE” S, ， 并 保证 S, 足够 大 且 S, 足够 小 ， 其 优化 目标 为 

min2A(4 XM XTA) + AA (10) — 
st. ATXHX' A-I argmaxJ (w) wr SW (14) 
其 中 : 47XZHX74=T ， 瑟 是 中 心 矩 阵 ， 工 是 单位 矩阵 ， M. 采用 拉 格 朗 日 乘 子 法 ， 即 可 求解 出 最 优 的 投影 向 量 : 


表示 由 源 域 样本 和 标签 ， 以 及 目标 上 样本 和 预测 伪 标 签 计算 w=S2(41w 一 轴 ) 。 在 实际 的 运动 想象 模式 识别 中 ， 本 文 将 学 
的 MAMD。 通 过 同样 计算 方法 ， 获 得 映射 矩阵 4 ， 用 于 变换 习 到 的 跨 域 特征 投影 至 最 优 方向 ， 并 采用 LDA 方法 进行 模 
源 域 和 目标 域 特征 。 式 识 别 。 
2.3 动态 分 布 自 适 应 3.8 运动 想象 EEG 信和 号 跨 域 特征 学 习 方法 流程 
JDA 方 法 将 边缘 分 布 和 概率 分 布 的 重要 性 视 作 同等 ， 然 本 文 提出 的 基于 MI-BCI 的 EEG 信号 跨 域 特征 学 习 方法 
而 ， 由 于 不 同 BCI 任 务 和 不 同人 群 ， 其 重要 性 并 不 始终 相等 。 流程 图 如 图 1 所 示 。 该 方法 步骤 如 下 : 
针对 JDA 的 不 足 之 处 ， 动 态 分 布 自 适应 方法 (BDA) 设 置 平衡 (1) 根据 不 同 个 体 域 的 EEG 样本 分 布 情况 ， 采 用 trial- 
妹子 uel 动态 调整 边缘 分 布 和 条 件 分 布 的 重要 性 031。 形 and-error 方法 选择 最 优 的 协 方差 距离 度量 方法 ， 随 后 ， 将 不 
式 地 ，BDA 的 优化 目标 为 同 个 体 域 的 EEG 信号 样本 对 齐 至 单位 阵 ， 使 得 个 体 域 的 
c > EEG 信号 样本 分 布 更 相似 ; 
mine [arx (o RE $m. Jra | utat (11) (2) 针对 每 个 域 对 齐 后 的 EEG 信号 样本 提取 CSP 特征 ， 
st. AXHX"A=1,0< ucl 在 此 基础 上 将 其 中 某 个 被 试 者 的 CSP 特征 作为 目标 域 ， 余 


其 中 : BDA 方 法 采用 A-distance 估计 。 需 要 指出 的 是 ，BDA ”下 所 有 被 试 者 的 CSP 特征 作为 源 域 ， 在 源 域 和 目标 域 间 执行 领 
方法 每 轮 和 迭代 后 ， 利 用 伪 标 签 和 A-distance 计算 新 的 4 值 。 域 自 适 应 方法 ， 使 得 CSP 特征 之 间 的 MD 最 小 化 ， 有 具体 来 看 ， 


同样 地 ， 获 得 的 映射 矩阵 用 于 变换 源 域 和 目标 域 特 征 。 J TCA, JDA, TJM 和 BDA 方法 中 选择 最 优 的 领域 自 适 应 方 
2.4 ”迁移 联合 匹配 法 ， 从 多 个 源 域 中 学 习 到 最 适合 目标 域 的 跨 域 特征 。 

无 论 是 JDA 还 是 BDA, 都 将 源 域 中 所 有 的 特征 整体 用 于 G) 采用 LDA 模型 用 于 跨 域 特征 分 类 ， 完 成 运动 想象 
计算 MMD ^ 从 而 获得 最 优 特征 映射 矩阵 4 。 实 际 上 ， 源 域 EEG 信号 的 模式 识别 。 


中 只 有 一 些 有 迁移 特征 意义 的 特征 ， 对 计算 MMD 能 够 起 到 HIA EE 

决定 作用 ， 而 另 一 些 是 存在 负 迁 移 的 特征 ， 则 会 对 MMD i 4 ”实验 与 结果 分 析 

算 造 成 反作用 。 因 此 ，TJM 方法 0 在 计算 放射 矩阵 时 ， 将 41 实验 初始 化 

源 域 的 放射 矩阵 添加 稀疏 正则 化 Aloa ， 而 保证 目标 域 的 放 为 了 验证 本 文 提出 的 跨 域 特 征 学 习 方 法 的 可 行 性 与 有 效 

射 矩 阵 为 二 范 数 不 变 lf 。 通 过 系数 正则 化 的 选择 ， 即 可 在 性 ， 本 文选 择 2008 BCI Competition IV 中 的 两 个 公开 的 MI- 

优化 过 程 中 从 源 域 中 选择 更 具 代 表 性 的 特征 ， 用 于 计算 BCI 数据 集 B637，MI 1 和 MI 2a。 该 两 个 数据 集 基 本 实验 范 

MMD ,从 而 提升 特征 迁移 的 性 能 S53。 因此 ， 在 TCA 基础 上 式 如 图 2 所 示 : 被 试 者 舒服 坐 在 椅子 上 ， 面 对 屏幕 开始 实验 。 

TJM 的 定义 可 写 为 t=0 时 ， 屏 幕 出 现 十 字 架 集中 注意 力 ， 伴 随 短 暂 提示 音 开始 
本 次 实验 ; t-2 时 ， 屏 幕 将 随机 出 现 指导 运动 想象 方向 的 箭 

(2) ， 头 ， 被 试 者 根据 箭头 指导 的 方向 进行 运动 想象 任务 : {=6 时 ， 


mintr (A"KMK TA)+A(lA; ka +A) 


Hc A RURAS I PORE Oe RH SS . MEE HS DU SE 
Ep: 4 为 正则 项 的 系数 。 由 于 上 述 优化 目标 具有 非 凸 的 “入 下 一 次 实验 。 
特征 ， 因 此 需要 通过 和 迭代 的 方式 求解 。 最 终 ， 由 4=[4.,4] 数据 集 MI I: 该 数据 集 采 集 了 7 名 健康 被 试 者 运动 想象 


构建 稀疏 映射 矩阵 ， 在 变换 源 域 和 目标 域 特征 的 同时 保证 稀 H EEG 数据 。 采 集 到 的 数据 包含 59 个 通道 ， 采 样 率 为 100 Hz. 
琉 性 ， 提 取 最 具 代 表 性 特征 。 本 文 实验 从 三 个 运动 想象 任务 中 提取 左手 和 右手 类 别 ， 分 别 标 
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录用 定稿 PAAA, $: 运动 想象 脑 电信 号 的 跨 域 特征 学 习 方 法 第 39 卷 第 8 期 
记 为 Classl 和 Class2， 每 个 用 户 的 每 个 类 别 有 100 个 样本 。 250Hz。 本 文 实验 从 四 个 运动 想象 任务 中 提取 左手 和 右手 类 


数据 集 MI 2a: 该 数据 集 采集 了 9 名 健康 被 试 者 运动 想 。 别 ， 分 别 标记 为 Classl 和 Class2， 每 个 用 户 的 每 个 类 别 有 72 
象 的 EEG 数据 。 采 集 到 的 数据 包含 22 个 通道 ， 采 样 率 为 。” 个 样本 。 


协 方差 矩阵 对 广 


I; 
Y 


EEG 脑 电信 号 
i 领域 自 适 应 方法 LDA 分 类 器 


协 方差 矩阵 对 章 


被 试 者 9 被 试 者 9 被 试 者 9 

图 1 跨 域 特征 学 习 流 程 图 
Fig.1 Cross-domain feature learning flow chart 

在 跨 域 特征 学 习 和 模式 识别 实验 中 ， 首 先 在 Matlab 上 均值 ， 分 类 效果 最 坏 的 是 调和 均值 。 与 MI 1 数据 集 的 分 析 

对 两 个 数据 集 进行 预 处 理 ， 将 运动 想象 提示 箭头 后 的 [0.5， ”过 程 相同 ， 第 2 个 被 试 者 在 调和 均值 下 的 分 类 准确 率 仅 有 

3.5] 秒 之 间 的 采样 点 数 作为 一 个 运动 想象 样本 。 因 此 ，MI 1 ”50.69%， 是 8 种 均值 下 表现 最 差 的 ， 因 此 导致 最 后 的 平均 结 


数据 集 每 个 样本 采样 时 间 点 数 为 300; MI 2a 数据 集 每 个 样 果 较 差 。 在 KL 散 度 均值 下 表现 最 差 的 分 类 准确 率 是 54.86%, 
本 采样 时 间 点 数 为 750。 比 其 他 均值 下 的 最 差 表 现 都 要 好 ， 单 个 最 差 值 对 KL 散 度 均 
实验 硬件 配置 为 Intel(R) Core(TM) i7-9700 CPU à  ” 值 的 影响 最 小 ， 因 此 KL 散 度 均值 的 效果 最 好 。 

3.00GHz，8GB 内 存 ， 运 行 64 位 Windows 10 系统 、 表 1 MI1 数据 集 的 不 同 均值 下 的 分 类 准确 率 
Matlab R2017b 以 及 Python 3.7.4。 协 方差 矩阵 中 心 对 齐 方法 Tab. 1 Classification accuracy under different mean values of MII dataset 
采用 代码 : https://github.com/hehe91/EA ， 领 域 自 适应 方法 mean riemann logeuclid logdet wasserstein ale harmonic kullback sym 
代码 则 使 用 代码 : https:// github.com/ jindongwang/ transfer 1 775 71.5 68 715 735 755 64 67 
learning / tree/master/ code。 针 对 JDA, BDA 和 TJM 的 迭代 2 771.5 74.5 76.5 735 80.5 7 715 7 
过 程 中 ， 都 将 迭代 次 数 设 置 为 10， 其 他 参数 则 针对 实际 数 3 755 715 705 74 70 73 71 70 
据 集 的 测试 结果 设置 最 佳 参数 。 本 文 取 所 有 迭代 结果 中 的 最 4 75 77.5 80 — 8L5 — 825 84 — 8 68 
佳 性 能 作为 最 终 的 分 类 结果 。 5 925 90 905 92 91.5 89 91.5 85.5 

6 755 77 765 7 7-5 75 75 74.5 

7 85 865 845 86 85 825 — 85 86.5 

Acc 79.70 7836 — 7807 7936 795 7936 7 74.93 


加 表 2 MI2a 数据 集 的 不 同 均值 下 的 分 类 准确 率 
soo o oea S o o 
短 


出 现 固定 交叉 


Tab.2 Classification accuracy under different mean values of MI 2a dataset 


mean riemann logeuclid logdet wasserstein ale harmonic kullback sym 


1 875 8472 8403 8405 8472 8333 8194 85.42 
0 2 3 E 5 6 7 iu 2 5625 52.78 5347 5417 5417 5347 50.69 58.33 

图 2 运行 想象 范式 流程 3 9861 97.92 9792 97.92 9792 9792 9653 97.92 

Fig.2 The flow chart of motor imagery paradigm 4 73.61 69.44 . 70.14 69.44 . 70.83 70.14 70.14 72.22 

4.2. 协 方差 矩阵 中 心 对 齐 方法 实验 结果 分 析 5 50 5417 5486 5625 5139 5486 56.94 54.86 
为 了 验证 不 同 均值 度量 计算 的 协 方差 矩阵 中 心 对 齐 方法 6 64.58 65.97 66.67 65.97 6597 6528 65.97 65.28 


在 不 同 数据 集 上 的 对 齐 效 果 ， 分 别 记 录 数 据 集 上 平均 准确 率 7 68.75 6944 6944 7014 6875 6875 70.4 68.06 
最 高 的 均值 度量 ， 以 及 每 个 被 试 者 域 最 佳 的 均值 度量 ， 在 后 8 89.58 90.97 9028 9028 9097 91.67 91.67 89.58 
续 实 验 中 分 别 验证 两 种 情况 下 的 最 佳 跨 域 特征 学 习 结 果 。 表 9 72.92 75.60 75.60 75.69 7639 77.78 7639 75.69 
1 和 表 2 是 分 别 给 出 了 数据 集 MII 和 MI 2a 的 实验 结果 。 其 Acc 73.53 73.46 — 73.601 73.77 7346 73.69 7338 74.15 
中 ， 每 行 表 示 使 用 8 种 不 同 均值 对 应 的 同一 被 试 者 的 分 类 准 4.3 ” 跨 域 特征 学 习 实 验 结果 分 析 
确 率 ， 每 列表 示 不 同 被 试 者 在 同一 均值 计算 的 协 方差 矩阵 中 十 对 不 同 域 对 齐 后 的 EEG 样本 ， 分 别提 取 CSP 特征 ， 
心 对 齐 下 的 分 类 准确 率 ， 最 后 一 行为 平均 结果 。 从 表 1 和 表 并 使 用 LDA 获得 分 类 识别 准确 率 。 由 于 样本 数量 较 少 ， 分 
2 的 实验 结果 中 可 以 看 出 : 别 对 每 一 个 数据 集 选 用 交叉 验证 法 ， 即 采用 留 一 法 依次 选 出 
针对 MI 1 数据 集 ， 分 类 效果 最 好 的 是 欧式 均值 ,效果 最 目标 域 样本 和 源 域 样本 ， 每 轮 实 验 依次 选择 一 个 被 试 者 的 样 
坏 的 是 KL 散 度 均值 。 由 于 调和 平均 值 对 单个 低 于 平均 值 的 本 集 为 目标 域 ， 余 下 的 样本 集 为 源 域 ， 每 一 个 被 试 者 都 会 被 
值 特别 敏感 ， 而 本 质 上 KL 散 度 均值 是 算数 和 调和 均值 之 间 ”选择 为 目标 域 ， 最 终 使 用 分 类 的 平均 正确 率 作为 实验 结果 的 
的 几何 平均 值 ， 在 KL 散 度 均值 下 第 1 个 被 试 者 和 第 4 个 被 度量 。 
试 者 的 分 类 结果 表现 都 很 差 ， 导 致 对 整个 平均 效果 变 差 。 在 4.3.1 实验 结果 
MI 2a 数据 集 下 的 测试 结果 中 ， 分 类 效果 最 好 的 是 KL 散 度 在 上 文 的 协 方差 矩阵 中 心 对 齐 方法 实验 结果 基础 上 ， 本 
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录用 定稿 韦 泓 好 ， 等 ; 运动 想象 脑 电 信号 的 跨 域 特征 学 习 方 法 第 39 卷 第 8 期 
文 分 别 在 两 个 数据 集 上 继续 测试 样本 对 齐 对 领域 自 适应 特征 表 6 MI 2a 数据 集 每 个 个 体 域 最 优 跨 域 特 学 习 准确 率 及 方法 
学 习 分 类 性 能 的 提升 效果 。 本 节 单 独 对 四 种 领域 自 适 应 方法 Tab. 6 Optimal accuracy of each subject in MI 2a dataset for the 
进行 实验 ， 以 及 在 样本 对 齐 后 提取 的 CSP 特征 基础 上 ， 分 cross-subject feature learning method 
别 进行 TCA, BDA, JDA, TJM 四 种 领域 自 适应 方法 和 协 ”subjects 均值 对 齐 准确 率 /% 特征 域 迁移 准确 率 /% ”最 优 参 数 。 运行 时 间 
方差 矩阵 中 心 对 齐 的 对 比 实验 。 为 了 更 加 明显 的 分 析 在 协 方 1 mean 87.5 JDA/TIM 9236 Dim=l,lambda=0.5  3.0s 
差 矩阵 中 心 对 齐 方法 下 四 种 领域 自 适应 方法 的 优 劣 ， 所 有 对 2  kullback sym — 5833 JIDA 60.42 Dim=4,lambda=0.1 2.96s 
比 均 在 欧式 均值 对 齐 获取 到 的 CSP 特征 上 进行 ， 表 3 和 表 4 3 mean 98.61 JDA 99.31 Dim=l,lambda=0.4 2.93s 
分 别 给 出 了 两 个 数据 集 在 不 同 特 征 学 习 方法 上 的 实验 结果 。 4 mean 73.61 JDA/TIM 77.78 Dim=2,lambda=0.1 2.78s 
Ae 3 和 表 4 中 的 结果 可 以 看 出 ， 在 平均 分 类 准确 率 方面 : 5 harmonic 56.94 JDA 59.03  Dim-3,lambda-0.9 3.1s 
(1) 在 四 种 领域 自 适应 特征 学 习 方法 的 对 比 下 ， 跨 域 特 6 logeuclid 66.67 TIM 72.22 Dim=3,lambda=0.1 2.74s 
征 学 习 的 准确 率 高 于 领域 自 适 应 特征 学 习 的 准确 率 ， 实 验 结 7  loget/harmonic 70.14 TJM 70.83 Dim-llambda-0.2 2.84s 
果 充 分 说 明 协 方差 矩阵 中 心 对 齐 方法 对 领域 自 适应 特征 学 习 8  ale/harmonic 91.67 JDA/TIM 95.14 Dim-3,lambda-0.2 2.75s 
方法 有 着 至 关 重 要 的 作用 以 及 对 降低 个 体 分 布 差异 的 有 效 性 。 9 ale 77.78 JDA 84.72 Dim-2,lambda-0.! 2.89s 
D 不 管 跨 域 特征 学 习 方 法 选择 哪 种 核 函数 ， 本 文 可 以 从 表 5 和 表 6 中 的 结果 中 可 以 看 出 : 
发 现在 MI 1 的 数据 集 上 TJM 的 方法 最 优 ， 在 MI 2a 的 数据 (1) 经 过 最 优 度量 方式 对 齐 后 的 样本 ， 其 CSP 特征 能 够 
RE JDA 的 方法 最 优 。 所 以 在 欧式 距离 样本 对 齐 基础 上 ， 在 最 优 域 自 适 应 方法 和 参数 下 获得 再 次 提升 ， 平 均 提 升幅 度 
TM 方法 更 适合 MU 数据 集 的 样本 分 布 ，JDA 方法 更 适合 ”为 2.43% 和 3.4096; 
MI2a 数据 集 的 样本 分 布 。 D 针对 不 同 数 据 集 的 不 同 个 体 域 ， 样 本 对 齐 的 度量 方 
(3) 与 缺少 样本 对 齐 的 领域 自 适应 特征 学 习 方 法 相 比 ， 式 对 最 终 识 别 准确 率 有 较 大 影响 ， 花 费 一 定时 间 代价 获取 最 
在 协 方差 矩阵 中 心 对 齐 方 法 的 预 处 理 条件 下 ， 四 种 领域 自 适 优 的 识别 准确 率 能 够 极 大 的 提升 模式 识别 性 能 ; 
Ps 应 特征 学 习 方 法 均 得 到 了 提升 。 其 中 ， 在 MI 1 Xm EX (3) 实验 结果 充分 说 明 ， 本 文 提出 的 方法 对 个 体 分 类 性 
er 用 EA+TJM 方法 提高 了 1435 9€, Æ MRa 数据 集 上 采用 能 的 提升 ， 面 向 个 体 域 数 据 分 布 差异 灵活 选择 合适 的 方法 ， 
C3 EA+JDA 方法 提高 了 3.09 %。 随 着 算法 的 发 展 在 实际 应 用 可 保证 高 性 能 和 重 棒 性 。 
© 表 3 MI 1 数据 集 在 不 同方 法 下 的 平均 准确 率 为 了 对 比分 析 本 文 提 出 算法 在 模式 识别 性 能 上 的 提升 ， 
e Tab.3 Average accuracy of MI 1 dataset at different methods 本 文 针 对 两 个 数据 集中 每 个 被 试 者 对 比如 下 种 情况 的 运动 想 
=f kena D 0TH AEE 跨 域 特征 学 习 (%) 象 识别 准确 率 : (1) 非 跨 域 的 识别 方法 (CSP+HLDA); (2) 不 经 
e EA TCA BDA JDA TJM EA-TCA EA*BDA EA+JDA EA+TIM 过 样本 对 齐 而 直接 使 用 领域 自 适 应 方法 (CSP+DA+LDA); (3) 
CN RBF 66.2966.14 67.5 68.86 8121 — 8100 82.43 8321 经 过 最 优 均 值 度量 样本 对 齐 后 的 识别 方法 (EA+CSP+LDA); 
N Linear 66.5 65.21 67.71 67.93 8121 80.35 82.28 82.93 (4) 本 文 提出 的 最 优 跨 域 特 学 习 方 法 (EA+CSP+DA+LDA)。 
K Primal 66.3666.71 67.36 68.21 81.57 79.64 83.00 83.14 图 3 分 别 给 出 了 两 个 数据 集 在 上 述 四 种 情况 下 的 运动 想象 识 
d Sam 66.7166.71 67.43 68.79 8107 — 8121 82.5 82.79 别 准 确 率 对 比 。 从 图 3 中 的 结果 可 以 看 出 ， 针 对 每 个 个 体 域 
2 表 4 MI 2a 数据 集 在 不 同方 法 下 的 平均 准确 率 的 最 优 跨 域 特征 学 习 方法 的 识别 准确 率 均 高 于 对 比 的 三 个 方 
x< Tab. 4 Average accuracy of MI 2a dataset at different methods 法 。 此 外 ， 最 优 度量 样本 对 齐 效果 基本 由 于 经 典 域 自 适 应 方 
a quu FEAR 。 领域 自 适应 (%) 跨 域 特征 学 习 (%) 法 。 对 比 非 跨 域 识别 方法 和 本 文 提出 的 最 优 跨 域 特征 学 习 方 
c EA TCA BDA JDA TJM EA*TCA EA+BDA EA+JDA EAHTJM 法 ， 运 动 想象 识别 准确 率 得 到 了 明显 提升 。 例 如 MI 1 数据 
rz RBF 72.0072.92 73.6174.07 73.99 — 7446 76.70 7647 集 的 subject4， 最 优 跨 域 特征 学 习 方法 比 非 跨 域 识别 提升 了 
O Linear as 72.5373.00 73.61 73.69 74.69 74.86 76.16 75.38 36.5%. 
Primal 72.2273.00 73.2373.69 74.31 74.46 75.62 76.00 同时 ， 如 图 4 所 示 ， 本 文 还 分 析 了 在 样本 对 齐 下 参数 
Sam 72.0072.92 74.1574.31 7438 | 7477 7608 76.62 dim 对 四 种 域 自 适应 特征 学 习 方 法 的 影响 。 从 图 4 的 结果 中 
从 两 个 数据 集 总 体 跨 域 特征 学 习 的 结果 可 以 看 出 ， 本 文 ” 可 以 看 出 ， 参 数 dim 的 选择 对 实验 影响 较 大 。 当 dim 较 小 


提出 方法 在 每 个 数据 集 上 均 取 得 了 比 单一 跨 域 特征 学 习 角度 。 时 ， 


所 能 承载 的 可 迁移 CSP Ti 因此 造成 ] 


识别 准 


任 较 少 ， 


更 好 的 运动 想象 模式 识别 性 能 。 为 了 探索 单个 个 体 域 特征 学 确 率 的 降低 。 随 着 dim 的 提升 ， 总 体 识别 准确 率先 上 升 后 下 
习 的 最 佳 方 法 ， 本 文 针 对 每 个 个 体 域 选择 最 优 的 样本 对 齐 度 降 ， 在 dim=4 时 获得 了 最 佳 的 分 类 性 能 。 
量 方法 和 CSP 特征 跨 域 特征 方法 ， 通 过 实验 获得 个 体 最 优 随后 ， 本 文联 合 跨 域 特征 学 习 的 实验 ， 选 择 rbf 核 函 数 ， 
的 运动 想象 模式 识别 结果 。 表 5 和 表 6 分 别 给 出 了 两 个 数据 “验证 在 它们 共同 作用 下 对 分 类 性 能 的 影响 ，MI1 和 MI2a Zt 
集 个 体 域 最 优 运动 想象 识别 性 能 的 方法 。 据 集 上 的 实验 结果 分 别 如 图 5 所 示 。 从 图 中 本 文 可 以 看 出 : 
表 5 MI1 数据 集 每 个 个 体 域 最 优 跨 域 特 学 习 准 确 率 及 方法 (1) 从 分 类 整体 情况 可 以 看 出 ，MI 1 数据 集 使 用 欧式 均 
Tab. 5 Optimal accuracy of each subject in MI 1 dataset for 值 作为 距离 度量 的 领域 自 适应 方法 效果 最 佳 ， 使 用 KL BUS 
the cross-subject feature learning method 均值 作为 距离 度量 的 效果 最 差 。MI 2a 数据 集 使 用 对 数 欧式 
subjects ”均值 对 齐 ”准确 率 /% 特征 域 迁 移 准确 率 /% ”最 优 参 数 。” ”运行 时 间 均值 作为 距离 度量 的 效果 比较 好 ， 对 比 协 方差 矩阵 中 心 对 齐 
1 mean 77.5 BDA 79 。 Dim=5,lambda=0.7 3.76s 方法 获得 的 分 类 准确 率 都 有 了 一 定 的 提升 ， 使 用 调和 均值 的 
2 wasserstein 805 JDA/TIM 845 Dim=2,lambda=0.7 3.60s 效果 整体 比较 差 。 
3 mean 75.5 TIM 79 Dim=l,lambda=0.5 3.63s (DMI 1 数据 集中 经 过 样本 对 齐 的 四 种 领域 自 适应 方 
4 ale 84  BDA/DA 865  Dim-6lambda-04 3.73s 法 中 效果 最 好 的 是 TIM. MI 2a 数据 集中 效果 最 好 的 是 
5 mean 92.5 TIM 95  Dim=l,lambda=0.5 3.81s JDA。 由 此 可 看 出 ， 采 用 欧式 均值 和 TIM 领域 自 适应 的 
6 ale 77.5 BDA 78.5 Dim=6,lambda=0.4 3.60s 特征 学 习 方 法 更 适合 MI 1 的 样本 分 布 ， 采 用 对 数 欧 式 均 
kullback sym 值 和 JDA 域 自 适 应 的 特征 学 习 方法 更 适用 于 MI 2a 数据 
7 86.5 TJM 88.5 Dim=2,lambda=0.1 3.51s 


riemann 


集 的 样本 分 布 。 
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图 3 四 种 不 同情 况 下 每 个 被 试 者 的 运动 想象 识别 准确 率 对 比 
Fig.3 Comparison of MI recognition accuracy of each subject under four different situations on MI 1 and MI 2a datasets 
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Fig.4 Comparison of accuracy of four domain adaptive methods in MI 1 and MI 2a data sets under different dim values 
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Fig.5 Comparison of four domain adaptive methods for MIl and MI 2a dataset under different mean calculations 
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在 实验 采用 的 四 种 领域 自 适应 方法 中 ，TCA 方法 为 无 变化， 本文 对 四 种 领域 自 适 应 特征 学 习 方 法 做 了 可 视 化 处 理 


监督 方式 且 不 需要 多 次 迭代 ， 而 BDA, JDA 和 TIM 三 种 算 


EE nf xxr EE o 


法 均 需 迭 代 获 得 伪 标 签 以 及 求 
出 了 不 同 迭 代 次 数 BDA、JDA 和 TIM 三 种 方法 在 不 同 核 函 
数 下 的 迭代 ; 
EE B AM 


从 


的 次 数 变 得 更 好 
RAKERA EN 10 次 ， 


迁 代 即 可 获得 与 最 终 迁 代 下 


接口 


M} 


场景 ， 本 文采 | 
量 的 迭代 次 数 即 可 者 


4.3.2 特征 可 视 化 
为 了 更 直观 的 观察 跨 域 特征 学 习 对 EEG 样本 和 特征 的 


E 


Kf 


pa 


解 系数 和 矩阵。 因此 ， 


图 6 给 


到 中 的 结果 可 以 看 出 ， 算 法 的 效 
。 此 外 ， 虽 然 迭 代 对 比 实 验 中 将 
实际 上 绝 大 多 数 情况 仅 需 5 次 


匹配 的 


"e 
均 


Hi EO PS RR 
的 识别 准确 率 。 


任 确 率 。 在 实际 脑 机 


学 习 方 法 ， 仅 需 少 


本 文选 择 


展示 MI 1 数据 集中 第 4 个 被 试 者 以 及 MI 2a 数据 


集中 第 8 个 被 试 者 的 EEG 信和 号 跨 域 特征 学 习 方 法 的 可 视 化 
结果 。 可 视 化 采用 t-SNE， 将 特征 投影 至 2 维 空间 并 展示 ， 
图 7 和 图 8 给 出 了 数据 可 视 化 结果 ， 图 中 的 红色 和 蓝 色 圆 点 
分 别 代 表 两 个 不 同类 别 。 

从 图 7 中 可 以 看 出 ， 进 行 样本 对 齐 之 后 ， 红 色 圆 点 和 蓝 


色 三 角 相 互 之 间 交 错 比较 多 ， 说 日 


经 过 样本 对 齐 的 领域 自 适应 方法 后 ， 红 色 


交错 也 少 了 ， 能 


有 分 类 


的 效果 不 是 很 好 ， 但 采用 


pan 


点 和 蓝 色 三 角 的 间隔 


LC 


上 较 明 显 的 区 分 出 样本 的 两 个 类 别 。 


从 图 8 本 文 可 以 看 出 ， 在 样本 对 齐 之 前 ， 红 色 圆 点 和 蓝 


色 


P 


En 


) 


bp 有 比较 好 的 分 类 效果 ， 但 经 过 样本 对 齐 的 预 处 理 后 


o 
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图 中 红色 


P, 


GP, Ri 运动 想象 脑 电 信号 的 跨 域 特征 学 习 方法 


到 点 和 蓝 色 三 角 的 交错 点 明显 少 了 ， 在 采 / 


协 方差 矩阵 中 心 对 齐 方 法 基础 上 ， 将 相同 被 试 者 不 同 


fi 


EZITA, £ 


度 的 提升 。 


不 同类 别 
不 同 被 试 者 EEG 信号 档 


IESE; 


HARIS EUER: 协 方差 矩阵 中 心 对 齐 方法 可 将 
昌 分 布 的 边缘 分 布 差异 减 小 ， 保 证 分 类 器 不 受 


84r 


82r 


74 


72 


色 圆 点 和 蓝 色 蓝 点 的 区 分 度 也 有 了 一 定 程 。” 类 别 信号 的 边缘 分 布 差异 减 小 ， 而 JDA 和 BDA 方法 则 不 但 
减 小 边缘 分 布 ， 更 减 小 了 标签 参与 的 条 件 分 布 差异 ， 因 此 色 
的 提升 平均 识别 准确 率 。TJM 方法 保证 了 特征 选择 
上 ， 可 以 被 选择 后 的 特征 
生 ， 且 让 不 同类 别 的 特征 分 离 的 差异 更 大 。 


在 可 视 化 的 结 


zs 


本 的 边缘 分 布 影响 。 此 外 ，TCA 7j 
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(a) MI 1 
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(b)MI 2a 
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Fig.6The diagram of accuracy variation under different iterations 


第 4 个 被 试 者 ,EA 前 


SAM EAS TCA. 


第 4 个 被 试 者 .EA+TIM 
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图 7 MI 1 数据 自 


B4 个 被 试 者 经 过 T-SNE 的 可 视 化 结 


Fig.7MI 1 dataset represents the visualization results of the fourth subject through T-SNE 
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图 8 MI2a 数据 集中 


第 8 个 被 试 者 经 过 T-SNE 的 可 视 化 结果 


Fig.8MI 2a dataset represents the visualization results of the fourth subject through T-SNE 
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4.4 RB E EE 
图 9 给 出 了 本 文 方法 的 运行 时 间 复 杂 度 对 比 结果 。 其 中 ， 
图 9(a) 展 示 了 不 同 均值 度量 平均 运行 时 间 。 从 该 柱状 图 的 结 
果 可 以 观察 到 ， 两 个 MI 数据 集中 在 欧式 均值 度量 下 的 运行 时 
间 最 短 ， 而 基于 AJD 的 对 数 欧 式 均值 运行 时 间 最 长 ， 前 者 快 
2-10 倍 。 此 外 ，MII1 数据 集 数据 量 更 大 ， 运行 时 间 更 入。 
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Fig.9Running time comparison of eight mean and 

four domain adaptive methods 

以 欧式 均值 度量 为 基础 ， 图 9 (b)(c) 给 出 了 不 同 域 自 适 

应 方法 运行 时 间 。 在 BDA, JDA 和 TIM 三 种 算法 中 迭代 次 
数 为 5 次 时 已 经 接近 最 优 解 ， 从 图 8(b)(c) 中 结果 可 以 看 出 ， 
TCA 的 运行 时 间 最 短 ， 相 对 而 言 BDA、JDA 和 TJM 的 运行 
时 间 相 差 不 大 ， 且 均 在 LS 秒 以 内 。 因 此 ， 所 提出 对 脑 电信 
号 识别 方法 的 响应 时 间 足 以 应 对 运动 想象 刺激 范式 的 流程 ， 
可 应 用 于 在 线 脑 机 接口 。 上 述 结果 表明 ， 最 优 跨 域 特征 学 习 
方法 的 平均 时 间 复 杂 度 在 可 接受 范围 内 ， 同 时 在 实际 的 MI- 
BCI 应 用 中 使 用 跨 域 特征 学 习 方 法 时 ， 需 根据 实际 情况 同时 
考虑 准确 率 和 时 间 复 杂 度 ， 选 择 适 合 场景 的 协 方差 中 心 均值 
度量 及 领域 自 适 应 方法 。 
目前 ， 现 有 脑 机 接口 系统 (如 外 骨骼 辅助 康复 ) 遭 过 新 被 
试 者 时 ， 往 往 采 集 大 量 样本 训练 模式 识别 模型 。 从 图 2 可 看 
出 ， 单 个 样本 至 少 需 要 4 秒 时 间 ， 且 通常 需要 200-300 个 样 
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本 才能 稳定 识别 。 采 用 本 文 提出 的 跨 域 特征 学 习 方法 ， 新 的 


被 试 者 仅 需 采集 少量 样本 (100 个 以 内 )， 即 可 通过 跨 域 特征 


学 习 使 用 已 有 使 用 者 的 脑 电信 号 样本 集 ， 且 随 着 使 用 者 增多 ， 可 


fs) 


新 的 被 试 者 所 采集 少量 样本 的 最 优选 择 ， 即 可 使 用 现 有 样本 集 获 


的 脑 电信 和 号 样本 集 也 越 多 ， 对 模式 识别 的 效果 也 越 好 。 此 外 ， 


得 更 高 的 模式 识别 性 能 ， 其 时 间 代 价 在 可 接受 范围 内 。 
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结束 语 
针对 运动 想象 个 体 域 EEG 信号 获取 代价 高 且 特 征 分 布 


差异 大 的 特性 ， 本 文 提出 了 一 种 运动 想象 脑 电 信号 的 跨 域 特 


征 学 习 方 法 ， 用 于 提升 跨 个 体 域 脑 电信 和 号 模式 识别 性 能 。 该 


方法 首先 选择 最 优 的 协 方差 距离 度量 方法 ， 对 齐 各 域 样本 集 
并 提取 CSP 特征 。 其 次 ， 在 CSP 特征 基础 上 ， 选 择 合适 的 
领域 自 适应 方法 ， 从 多 源 域 中 学 习 出 最 适合 目标 域 的 跨 域 特 


征 ， 


所 提出 的 方法 明显 可 提升 模式 识别 准确 率 。 跨 域 特征 学 习 能 
够 极 大 的 提升 BCI 应 用 中 的 样本 数量 和 质量 ， 提 升 控制 外 部 
设备 的 性 能 ， 然 而 现 有 方法 存在 一 个 缺点 : 采用 离线 方式 进 


最 后 采用 LDA 分 类 完成 模式 识别 任务 。 实 验 结果 表明 ， 


行 。 


il 


在 实际 应 用 中 ， 源 域 样本 或 许 能 够 长 期 积攒 ， 但 是 目标 


域 样本 不 断 产生 ， 无 法 直接 计算 协 方差 矩阵 中 心 以 及 CSP 


FE。 今后 的 工作 集中 在 : 探索 从 代 式 协 方差 计算 方法 ， 设 


计 增 量 迭 代 式 的 跨 域 特征 学 习 方 法 ， 并 将 迭代 式 跨 域 特征 学 


JM 


至 在 线 MI-BCI 中 。 
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